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Abstract 

We investigate the effects of quenched randomness on the universal proper- 
ties of a two-temperature lattice gas. The disorder modifies the dynamical 
transition rates of the system in an anisotropic fashion, giving rise to a new 
fixed point. We determine the associated scaling form of the structure factor, 
quoting critical exponents to two-loop order in an expansion around the up- 
per critical dimension d c = 7. The close relationship with another quenched 

disorder fixed point, discovered recently in this model, is discussed. 
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The statistical mechanics of non-equilibrium phenomena remains poorly understood, in 
spite of the ubiquity of such systems in nature. As a first inroad towards the full complexity of 
such systems, simple model systems driven into non-equilibrium steady states have attracted 
much interest over the past decade [[]. To prevent the systems from settling into the usual 
Boltzmann equilibrium steady state, a uniform current (of particles or energy) is maintained, 
e.g., by imposing a uniform drive, or by coupling the system to two reservoirs at different 
temperature or chemical potential. The prototype model is an Ising lattice gas, in the 
presence of a uniform "electric" field acting on the "charged" particles [@]. Motivated by the 
physics of fast ionic conductors f3[, the usual Ising energetics is modified by the field which 
enhances (suppresses) the rate for nearest-neighbor particle-hole exchanges along (against) 
a specific lattice axis. The total number of particles is conserved, and periodic boundary 
conditions are imposed to ensure a non- vanishing global current. The resulting steady state 
distribution is non-Hamiltonian and violates detailed balance. In spite of its simplicity, this 
model, and many of its variants, exhibit surprising features like anomalous correlations at 
all temperatures and different types of non-equilibrium universal behavior JT|. 

In the framework of renormalization group (RG) studies of these models, the universal- 
ity class of order- disorder transitions is determined by the symmetries of the drive and the 
local order parameter alone, provided the dynamics is short-ranged in space and time. For 
equilibrium systems it is well known that spatial or temporal nonlocalities in the Hamilto- 
nian, such as sufficiently long-ranged interactions or quenched impurities, can have profound 
consequences for critical behavior |4]|| . On the other hand, non-local terms appearing only 
in the dynamics but not in the Hamiltonian have no effect on static universal properties 
provided detailed balance holds. In contrast, if detailed balance is broken, the steady state 
is determined by both the local energetics and the dynamical rules, and nonlocalities in 
either are expected to play an important role. 

So far, only a few studies have addressed this issue. The effects of long-ranged random 
spin exchanges, added to the non-conserved Glauber dynamics of the standard Ising model, 
have been investigated in some detail 0. The resulting steady state, while Hamiltonian, is 
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characterized by effective long-range interactions between spins. Other studies have focused 
on the effect of quenched disorder in the uniformly driven Ising lattice gas |7||1- Here, 
due to the presence of a global current, a subtle complication arises when attempting to 
compare field theoretic results J7| with Monte Carlo data ||. The simulations are performed 
on a finite lattice with periodic boundary conditions, so that particles typically see a given 
disorder configuration repeatedly in the course of thermal averaging, thus inducing effective 
correlations in the disorder average. In contrast, the analytic work rests on the assumption 
of uncorrelated randomness. 

Motivated by an interest in deepening our understanding of quenched disorder in driven 
lattice gases, while recognizing the need to establish a clear correspondence between Monte 
Carlo data and field theory results, a non-equilibrium lattice gas without global current was 
studied in ||. In the absence of quenched disorder (henceforth called the pure case), the 
field theory for this system captures the long-wavelength, long-time behavior of two different 
microscopic models: a randomly driven Ising model [TIJ and a two-temperature lattice gas 



Adding a quenched contribution to the drive in the randomly driven system, the authors 
showed that the field theory of the pure model is significantly modified, giving rise to a novel 
fixed point. The associated universal scaling behavior was investigated to two- loop order 
and found to be distinct from both the pure and the non-driven case. 

In this Letter, we continue the analysis of quenched disorder in this system. Modelling 
the disorder as a quenched random bias acting along longitudinal or transverse bonds, we 
show that these two scenarios give rise to two distinct fixed points. The former of these, 
subsequently referred to as the "longitudinal" fixed point, was considered in |J. Here, we 
investigate the other ("transverse") fixed point. We begin by summarizing the relevant 
properties of the pure system, followed by the introduction of quenched randomness. Next, 
we discuss the resulting effective Langevin equations for the longitudinal and transverse fixed 
points. The universal scaling behavior associated with the latter is then analyzed, using 
power counting techniques confirmed by an explicit e-expansion about the upper critical 
dimension of the theory. Explicit results for critical exponents are computed to two-loop 



order. We conclude with some comments. 

In the absence of frozen disorder, our field theory describes the coarse-grained properties 
of two different microscopic models. The first, a randomly driven system ||10|| , corresponds 
to a (half-filled) lattice gas of particles (up-spins) and holes (down-spins) on a fully periodic 
lattice, interacting through an attractive nearest-neighbor Ising coupling and coupled to a 
heat bath. As in the uniformly driven system, particle-hole exchanges along a specific lattice 
axis are biased by a driving field. Here, however, the amplitude E of the drive switches sign 
and magnitude randomly in space and time, i.e., a new value of E is culled from a symmetric 
distribution, p(E), for each pair exchange along the field direction. Exchanges along the 
transverse direction are controlled by energetics alone. The symmetry of p(E) ensures that 
there is no global current; otherwise its detailed form is not important for critical behavior. 
In its other version, a two-temperature model, we again consider a half-filled Ising lattice gas 
with spin-exchange dynamics. The drive, however, is now replaced by a second heat bath, 
such that spin pairs forming a "longitudinal" bond exchange according to a temperature 
Tji while "transverse" bonds are controlled by a different temperature T± < T» \TL\. Both 



microscopic models obey the same set of global symmetries, namely: translation invariance 
along the longitudinal and the transverse directions, as well as rotation invariance in the (d — 
l)-dimensional transverse subspace. Moreover, both local order parameters are scalar and 
exhibit the usual up-down symmetry of the Ising model. The other common characteristic 
is that, as the (transverse) temperature is lowered, long- wavelength fluctuations with purely 
transverse wave vector become soft first, indicating that criticality sets in when the transverse 
diffusion coefficient vanishes while the longitudinal one is still positive. A Langevin equation 
for the local magnetization (j)(x,t), consistent with these features, is easily postulated, or 
derived from its counterpart for the uniformly driven system by averaging with weight p(E) 

= ( Tj _ - v 2 )V 2 + th<9 2 + |V 2 3 - d n - V 3± . (1) 
Here d (V) denotes the gradients in the parallel (transverse) directions. j\\ and j± are 



Gaussian white noise terms which are distributed according to exp{— J d d xdt [j\ 2 /4n\\ + 
j± 2 /4n±\} and model the effect of thermal noise after coarse-graining. The parallel and 
transverse diffusion coefficients are denoted by rii and t±. Criticality is associated with 
t± — > while T|| > 0. Since r\\ does not vanish, T\\d 2 4> is the only relevant term involving 
d 2 . Otherwise, all coefficients here are functions of the microscopic control parameters but 



their explicit form is not required. Standard field theory methods [T2| can be employed 



to compute the universal scaling behavior of the theory HlQfl , in an expansion around the 



upper critical dimension d c = 3. High precision Monte Carlo data for the two-temperature 



model [13| are consistent with the field theory predictions, confirming that the two models 
do indeed belong to the same universality class. 

Next, we introduce frozen disorder into the microscopic dynamics of the system, either 
in its random drive or two-temperature version, by adding a quenched random bias A a (x) 
to the exchange rate of a particle-hole pair located on a bond in the a-direction, a = ||, _L. 
This type of randomness is easily implemented in a Monte Carlo simulation, since it can 
be considered as a time-independent local drive resident on each bond. In the Langevin 
equation (1), it leads to an extra term, V • J, on the right hand side, with the Ohmic current 
J a given by J a = c(0)A a (x), in analogy to the uniformly driven case QTJ] . Here, V is the 
full c?-dimensional gradient, and c(<f>) oc 1 + O(0 2 ) is the (density-dependent) conductivity. 
This form of J a is also easily derived from a microscopic hopping model, upon taking a 
naive continuum limit and neglecting higher order derivatives. These, as well as the O(0 2 ) 
correction to the conductivity, can be shown to be irrelevant in the RG sense. The A a (x) are 
taken to be Gaussian distributed with zero mean, < A a (x) >= so that no global current is 
induced, and < A a (x) Ap(x') >= 2a a 8 a p8(x—x'). The positive parameters a± and a\\ control 
the strength of the randomness. It is now apparent that the leading modification of (1), apart 
from trivial redefinitions of its coefficients, is a quenched random current, — dA\\ — VA_l, 
on the right hand side. To proceed, it is most convenient to introduce a Martin-Siggia-Rose 
field 4>{x,t) [|14|| , and to recast the Langevin equation as a dynamic functional [13]. In this 



form, the average over the A a (x) is easily performed, leading to an effective functional 



J [0, $\ = J d d x J dt\4> {a~ x - (r ± - V 2 )V 2 + r||9 2 + |v 2 3 J (2) 
+ Jd d xJ dtX J dt'X$(x, t) {{n\\d 2 + n±V 2 )X~ 1 8(t - t') + (a {l d 2 + a ± V 2 )} 0(x, t') . 

We emphasize that both the frozen disorder and the thermal noise generate operators 
quadratic in 0. In contrast to the thermal noise, however, the quenched random contribu- 
tion to the Langevin current leads to an operator which is nonlocal in time. Correlation and 
response functions can now be computed as functional averages with weight exp — J [0, <j)]. 

The next step in the RG treatment of (2) is a dimensional analysis. Since t± vanishes 
first, parallel momenta scale as k\\ oc fi 2 , on a characteristic scale \i oc k±. Thus, strongly 
anisotropic scaling QXJ] is observed already at the tree level. As a consequence, the "trans- 
verse" coefficient n± is more relevant than its "longitudinal" counterpart, n», which will 
therefore be dropped. Similarly, o~± , if non-vanishing, is more relevant than mi . A signif- 
icant advantage of our particular implementation of disorder, as a bias acting on parallel 
or transverse bonds, now becomes apparent: it allows us to tune a± or o~\\ separately to 0. 
Here, we will choose a± ^ 0, in contrast to Ref. 0. We will return to the consequences 
of either choice in the conclusions. Proceeding with the dimensional analysis, we can now 
neglect cry as well. Requiring that the remaining terms in (2) be dimensionless, we find 
Xt oc yU -4 , so that the thermal noise becomes irrelevant and therefore negligeable compared 
to the quadratic term associated with the frozen disorder. Furthermore, choosing the scale 
of 4> such that crj_ oc we find oc /i^ -5 ^ 2 and <fi oc fi^ d+7 ^ 2 , implying that g oc jj7~ d . The 
upper critical dimension follows as d c = 7, in contrast to the pure system where d c = 3. We 
also recall the result d c = 5 for a disordered system with o~± = [[J. Thus, we may expect 
novel universal behavior for our theory. We should remark that such dimensional shifts have 
also been found in the uniformly driven system with frozen disorder 0. 

We now consider the critical theory (tj_ = 0) with insertions of Ar_|_0V 2 0, to two-loop 
order in e = d c — d. A rescaling of allows us to set o~± = 1. The anisotropic forms 
of the bare correlation and response propagators are easily obtained from the Gaussian 
parts of (2): with A(k) = k]_(k]_ + t±) + T\\k 2 , we find G D (k,uj) = [iu + XA(k)}^ 1 for the 
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response propagator, and S (k,u) = An\ 2 8(ui) \iu + AA(fc)|~ 2 for the correlator. The latter, 
in particular, remains proportional to 5 (id) at all orders, due to the presence of the time- 
delocalized noise in (2). A simple scale invariance of the theory, under xu — ► 7x11 at constant 
x± with arbitrary scale factor 7, identifies the effective dimensionless expansion parameter 

— 1/2 

of the theory as u = ^ e T\\ 9- Denoting one-particle irreducible vertex functions with n 
(h) external 0- (</>-) legs and m insertions by , we find that only Tf} , T^, and 
are primitively divergent. Thus, only the quantities 0, u, and t± must be renormalized to 
render the theory finite. The coupling constant u flows towards an infrared stable fixed 
point u* , while the renormalizations of and t± determine the only two independent critical 
exponents, 77 and v. Following standard methods ||12|| , the renormalization group equations 
at the fixed point, combined with dimensional analysis and scale invariance, predict the full 
asymptotic scaling behavior of the vertex functions r~™\ Instead of quoting the complete 
expressions here, we specialize to the steady-state structure factor, which is easily measured 
in simulations: 

S(k, t; r ± ) = ^S(k\\/^, hjfi, *//*; rj^) (3) 

Some comments are in order here. First, the prefactor /i _6+r? may appear somewhat surpris- 
ing, compared to the more familiar fi~ 2+v . However, its origin is easily traced back to the 
tree level result for the structure factor, S a (k, id), where the S(u) is responsible for changing 
the familiar power of 2 to the more unusual 6. Second, the strong anisotropy exponent A 
|U and the dynamic exponent z are related to the two independent exponents rj and v via 
scaling laws: A = 1 — 77/4 and z = 4 — 77/2. The latter differs from the familiar z = 4 — 77 
(for the conserved Ising model), again by virtue of the quenched disorder which modifies 
the usual relationships between the Wilson functions for 0, 0, and A. The details of the 
two-loop calculation are too technical to be reported here, but the stability of the fixed 
point is established and explicit results for the critical exponents can be computed. We find 
77 = ^ffyC 2 + 0(e 3 ) and v = \ + ^ + ^ (^fp + In § + 6 In 4) + O (e 3 ), which shows clearly 
that our theory falls into a new universality class. 
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Other scaling laws determine the remaining anisotropic exponents 0, such as the two 
different correlation length indices u± — v and i/m = z/(l + A), or the two different dynamic 
exponents z± = z and z n — z/(l + A). Due to the anisotropy, care must also be taken when 
defining exponents that characterize the power-law decays of the critical structure factor, in 
real and momentum space. Simple scaling laws relate these indices to i] and A Finally, 
the order parameter exponent (3 = \v{d — 6 + r\ + A) follows from the scaling form of the 
equation of state ]12] . 

In summary, we have analyzed the critical behavior of a two-temperature stochastic 
lattice gas with conserved dynamics, in the presence of a quenched random bias acting 
on bonds in both the longitudinal and transverse subspace. The bias generates a time- 
delocalized Langevin current which modifies the usual power counting so dramatically that 
the upper critical dimension is shifted to d c = 7. The RG analysis yields the asymptotic 
scaling forms of correlation and response functions; and critical exponents, associated with 
a new universality class, have been computed to 0(e 2 ) in an expansion about d c . If, in 
contrast, the frozen disorder affects longitudinal bonds only, then a± = in Eqn (2). In 
this case, the dimensional analysis has to be revisited, resulting in d c = 5. Moreover, an 
additional relevant operator — k J d d x J dt J dt'\ 2 (j)(x,t)'V 4l <f)(x,t') is generated under the RG 
and must be included in the dynamic functional. Considering the interplay between the two 
couplings (Tii and k, the former is found to be irrelevant while the latter gives rise to a new 
( "longitudinal" ) fixed point which is infrared stable near d = 5 0. It is noteworthy || that 
this fixed point also captures the effect of quenched disorder resident in the nearest-neighbor 
hopping rates, modelled, e.g., by a Gaussian distribution of potential barriers between sites. 
Eqn (2), including the operators associated with cry, a± and k thus provides a unified field 
theory for several different types of quenched disorder. Work is in progress to clarify the 
intriguing connections between different microscopic implementations of such disorder and 
their surprising degree of universality (16| . 
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